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We propose a new method to determine at what transverse momenta particle production in 
high-energy pp collisions is governed by hard parton-parton processes. Using information on the 
transverse spatial distribution of partons obtained from hard exclusive processes in ep/'yp scattering, 
we evaluate the impact parameter distribution of pp collisions with a hard parton-parton process 
as a function of pr of the produced parton (jet). We find that the average pp impact parameters 
in such events depend very weakly on pr in the range 2 < pr < few 100 GeV, while they are 
much smaller than those in minimum-bias inelastic collisions. The impact parameters in turn 
govern the observable transverse multiplicity in such events (in the direction perpendicular to the 
trigger particle or jet). Measuring the transverse multiplicity as a function of pr thus provides an 
effective tool for determining the minimum pr for which a given trigger particle originates from a 
hard parton-parton process. Additional tests of the proposed geometric correlations are possible by 
measuring the dependence on the trigger rapidity. Various strategies for implementing this method 
are outlined. 
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I. INTRODUCTION 

The first experimental results from LHC once again 
raise the question at what transverse momenta particle 
production in pp collisions is dominated by hard parton- 
parton interactions. A quantitative understanding of the 
relevant mechanisms is important not only for future 
studies of QCD phenomena, but also for controlling the 
strong interaction background in new particle searches. 
The challenge lies in the fact that the growth of the av- 
erage multiplicities makes it very difficult to observe jets 
with moderate pt, while at the same time the properties 
of non-perturbative semi-hard dynamics and its ability 
to produce particles with px ~ few GeV are not well 
understood. 

In an earlier article we demonstrated that the nu- 
cleon's transverse partonic structure plays an essential 
role in the theoretical analysis of pp collisions with hard 
processes. Experiments in hard exclusive electroproduc- 
tion of vector mesons "f*p — > V + p and photoproduction 
of heavy quarkonia jp — > J ftp + p have shown that the 
gluons with 10~ 4 < x < 10 _1 are localized at small trans- 
verse distances of 0.4 — 0.5 fm (median, depending on x 
and Q 2 ), much smaller than the characteristic range of 
soft interactions at high energies, see Fig. QJ,. Qualita- 
tively, this is explained by the fact that Gribov diffusion 
in the partonic wave function, which causes the range 
of soft interactions to grow with energy Q , is suppressed 
for highly virtual constituents. In pp scattering this two- 
scale picture implies that hard processes mostly occur in 
central collisions, where the areas occupied by partons in 
the relevant x-range overlap. Peripheral collisions con- 
stitute the dominant part of the overall inelastic cross 
section without contributing much to inclusive jet pro- 



duction, see Fig. [lb. A trigger on a hard process thus, on 
average, selects central pp collisions [![. Numerical stud- 
ies show that at a center-of-mass energy y/s — 14TeV 
a dijet trigger on pr ~ few 10 GeV reduces the median 
pp impact parameter b by a factor of ~ 2 compared to 
minimum-bias inelastic collisions; the reduction is nearly 
as strong at the current LHC energy of 7 TeV (see below) . 

Here we point out that these insights into the trans- 
verse geometry of pp collisions can be used to address the 
question at what transverse momenta particle production 
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FIG. 1. (a) The two-scale picture of transverse nucleon struc- 
ture at high energies (transverse view), (b) Its implication 
for pp collisions. Peripheral collisions constitute the domi- 
nant part of the overall inelastic cross section. Hard processes 
happen predominantly in central collisions, where the areas 
occupied by large-a: partons overlap. 
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is governed by hard parton-parton processes. The key 
observation is that the transverse multiplicity, measured 
in the direction perpendicular to the transverse momen- 
tum of the trigger particle (or jet), is correlated with the 
average impact parameter in the underlying pp event. If 
the trigger particle originates from a hard parton-parton 
process, the underlying pp collision is central, and the 
average impact parameter depends only weakly on px- 
The transverse multiplicity thus should be practically in- 
dependent of pt, and substantially larger than that in 
minimum bias events. If, however, the trigger particle is 
produced by soft interactions, the transverse multiplicity 
should be substantially smaller and reflect the average 
multiplicity in minimum-bias inelastic collisions. In this 
sense, the transverse multiplicity can serve as a diag- 
nostic of the dynamics in the production of the trigger 
particle at given pt- First theoretical suggestions for 
studies of the transverse multiplicity were put forward 
in Ref . 3] . Experimental investigation of the correlation 
between the jet production and the structure of the un- 
derlying event were pioneered by the CDF experiment 
at the Tevatron The first data on underlying event 
structure in collisions with hard processes at the LHC 
were recently reported by ATLAS and CMS HQ. 

Additional tests of the geometric correlations described 
here become possible with measurements of the depen- 
dence of the transverse multiplicity on the rapidity of the 
trigger jets. In particular, we predict that in the rapidity 
region not affected by the fragmentation of the trigger 
jets the enhancement of the multiplicity will persist and 
be isotropic in transverse space. In this way one could 
verify the universality of particle production in the cen- 
tral pp collisions selected by hard processes. 

This article is organized as follows. In Sec. [II] we sum- 
marize our knowledge of the nucleon's transverse partonic 
structure and update our paramctrization of the trans- 
verse gluonic size as a function of x. In Sec. IIIII we use 
this information to study the impact parameter distri- 
butions of pp events with hard processes in dependence 
on the transverse momentum px of the trigger particle 
(jet), for the kinematics currently covered by LHC, and 
find that the median b weakly depends on the trigger pt 
and is substantially smaller than that in minimum-bias 
inelastic collisions. In Sec. IIVI we discuss the connec- 
tion between centrality and the transverse multiplicity, 
and how measurement of the latter provides an effective 
means of quantifying at what pr particle production is 
dominated by hard QCD processes. In Sec. [V] we con- 
sider the dependence of the multiplicity on the rapidity 
of the trigger, and how it can be used for additional tests 
of the dominance of small impact parameters in pp colli- 
sions with hard processes. In Sec. IVII we present several 
suggestions for further analysis of the pp event structure 
data. A summary and discussion of our results are given 
in Sec. rVTH 

Our analysis relies on information on the nucleon's 
transverse partonic structure obtained from hard exclu- 
sive processes in ep/^yp scattering. Extending our previ- 



ous study we present here an updated parametriza- 
tion of the transverse distribution of gluons, which takes 
into account the more recent HERA data @, ©] and per- 
mits realistic uncertainty estimates. The numerical re- 
sults are nevertheless close to those obtained in our pre- 
vious study. 

Current Monte Carlo (MC) generators for pp events 
usually do not take into account the available experi- 
mental information on transverse nucleon structure and 
treat the distribution of gluons over transverse position 
as a free function. The typical setting for PYTHIA 
and HERWIG correspond to a transverse area occu- 
pied by gluons which is a factor ~ 2 smaller than what is 
indicated by the HERA data (see below). In the analysis 
of experimental data, the shape of the transverse gluon 
distribution is usually treated as one of the tuning pa- 
rameters, see e.g. Refs. 0, H|. While we do not directly 
address these technical issues here, our results certainly 
have implications for the design of future MC generators 
for pp events at LHC. 



II. TRANSVERSE PARTONIC STRUCTURE OF 
NUCLEON 

Information on the transverse spatial distribution of 
gluons in the nucleon comes from the study of hard exclu- 
sive processes such as electroproduction of vector mesons, 
7*p — > V + p, or the photoproduction of heavy quarko- 
nia, "fp — t J /ip + p. Thanks to a QCD factorization 
theorem [l2| . the amplitude of these processes in the 
leading-twist approximation can be expressed in terms 
of the gluon generalized parton distribution (or GPD), 
which parametrizes the matrix element for the emission 
and absorption of a gluon by the target. Of particular in- 
terest is its i-dependence in the "diagonal" case of equal 
momentum fraction of the emitted and absorbed gluons. 
It is described by the normalized two-gluon form fac- 
tor F g (x, t\Q 2 ), where x is the gluon momentum fraction 
and t = — Aj^ the transverse momentum transfer to the 
target. This function can be regarded as the transverse 
form factor of gluons with longitudinal momentum frac- 
tion x in the nucleon. Its Fourier transform describes the 
transverse spatial distribution of gluons with given x, 

F g (x,p\Q 2 ) ee ||^ e KAxp) Fg{x ,t = -Ai\Q% 

w 

where p = \p\ measures the distance from the transverse 
center-of-momentum of the nucleon, and the distribu- 
tion is normalized such that / d 2 p F g (x , p\Q 2 ) = 1. 

Experiments in hard exclusive processes actually probe 
the gluon GPD in the non-diagonal case (different mo- 
mentum fraction of the emitted and absorbed gluon) , be- 
cause of the longitudinal momentum transfer required by 
kinematics. At x < 10~ 2 , QCD evolution allows one 
to relate the non-diagonal gluon GPD to the diagonal 
one at the input scale, and the diagonal two-gluon form 
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factor can be directly inferred from the f-dependence of 
the measured cross sections. At larger x, the relation 
between the non-diagonal and diagonal GPDs generally 
becomes model-dependent, but useful information can 
still be extracted using GPD paramctrizations. 

The i-dependence of the measured differential cross 
sections of exclusive processes at \t\ < 1 GeV 2 is com- 
monly described either by an exponential, or by a 
dipole form inspired by analogy with the nucleon elas- 
tic form factors. Correspondingly, we consider here two 
parametrizations of the two-gluon form factor: 



F g (x,t\Q 2 ) = 



exp(B g t/2) 

(1 - t/mly 



(2) 



where the parameters B g and m g are functions of x and 
Q 2 . The two parametrizations give very similar results 
if the functions are matched at \t\ = 0.5 GeV 2 , where 
they are best constrained by present data (see Fig. 3 of 
Ref. 13]); this corresponds to 



B g = 3.24/m 2 . 



(3) 



We use both parametrizations in our studies below; the 
difference between the results serves as an estimate of the 
uncertainty due to our lack of precise knowledge of the 
shape. The corresponding spatial distributions of gluons 
in the transverse plane, Eq. (TJJ), are given by 



F g (x,p\Q 2 



(27TB,)- 1 exp[-p 2 /(2 J B fl )] ) 



K/(2tt)] (m g p/2) K^mgp), 



(4) 



where K\ denotes the modified Bessel function. 

Most of the experimental information on the nucleon's 
two-gluon form factor comes from J ftp photoproduc- 
tion, which probes the gluon GPD at an effective scale 
Q 2 ~ 3 GeV 2 , determined by the average transverse size 
of the cc pair during its interaction with the target, and 
momentum fractions of the order x ~ M 2 ,^/W 2 [lij ]. 
When extracting the two-gluon form factor from the 
slope of the J/tp differential cross section, a correction 
is made for the effect of the finite J/tp size on the ob- 
served i-dependence, 



B 



7/V 



= B„ 



AS, 



(5) 



where AB s» 0.3 GeV from a dipole model estimate [14j . 

The data from the FNAL E401/E458 broadband beam 
experiment at (x) — 0.05 [l5| . in which the recoiling pro- 
ton was detected, are described by an exponential two- 
gluon form factor with B g = 3.0 GeV -2 (see Fig. [2]), al- 
beit with large errors, or a corresponding dipole form 
factor with m 2 — 1.1 GeV 2 . Comparison with the mass 
parameter in the dipole parametrization of the nucleon's 
electromagnetic (Dirac) form factor, m 2 m = 0.7 GeV 2 , 
indicates that at these values of x the average transverse 
gluonic radius squared is only ~ 0.6 of the transverse 
electromagnetic radius squared. 
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FIG. 2. The exponential t-slope, Bjm, of the differential 
cross section of exclusive J/tp photoproduction measured in 
the FNAL E401/E458 0, HERA HI [<|, and ZEUS |rJ 
experiments, as a function of x = Mj^/W 2 . (In the HI 
and ZEUS results the quoted statistical and systematic un- 
certainties were added linearly.) The dashed lines represent 
the published two-dimensional fits to the HI and ZEUS data 
[|,[fj]. The parameter B g in the exponential two-gluon form 
factor Eq. ([2]) is related to the measured J/ip slope by Eq. 1(5). 
Our parametrization Eqs. ©-I© is shown by the solid line. 



The HERA data at x < 10~ 2 @, [l6| indicate that the 
transverse gluonic size of the nucleon increases with de- 
creasing x (see Fig. [5]). In the region x ~ fewlO -2 this 
effect is partly due to the contribution of the nucleon's 
pion cloud, which is strongly suppressed at x > 0.1 and 
fully developed only for x < 0.01 [l7[- Another con- 
tribution may arise from Gribov diffusion, which is sup- 
pressed by the hard scale Q 2 but still not negligible. Over 
the range covered by HERA, the rate of increase of the 
gluonic size can be parametrized by an effective Regge 
slope, a' g , which is smaller than that for corresponding 

soft processes, a' soft = 0.25 GeV -2 . Averaging the fits to 
the HERA HI and ZEUS data @, OH, and accounting for 
the correction Eq. ((5)), we parametrize the x-dependence 
of the gluonic exponential slope as (here Q 2 = 3 GeV 2 ) 



B g (x) = B g0 + 2a' g \n(x /x), 
x = 0.0012, 
B a o = 4.1 (±S:1) GeV" 2 , 
a\ =0.140 GeY" 2 . 



(6) 
(7) 
(8) 
(9) 



The uncertainties in parentheses represent a rough es- 
timate based on the range of values spanned by the HI 
and ZEUS fits, with statistical and systematic uncertain- 
ties added linearly. One sees from Fig. [5] that the fit to 
the HERA data consistently extrapolates to the FNAL 
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data point. The corresponding dipole parametrization 
obtained via Eq. © is close to the one used in our pre- 
vious study [![. 

The transverse spatial distribution of partons also 
changes with the resolution scale, Q 2 , as a result of 
DGLAP evolution. Generally, the partons observed at 
a given momentum fraction x and scale Q 2 are decay 
products of partons with x' > x which existed at a lower 
scale, Qq. In the leading-twist approximation the de- 
cay happens locally in transverse space. As a result, the 
transverse size observed at fixed x shrinks with increas- 
ing Q 2 , because the decaying partons at the lower scale 
had larger momentum fractions and were localized in a 
smaller transverse area. In order to calculate the change 
of the transverse spatial distribution of gluons with Q 2 
one would need to know the spatial distributions of both 
gluons and singlet quarks at the input scale for all x' > x, 
where they are only poorly constrained by present data. 
Numerical studies based on a simple parametrization [TJ 
suggest that the evolution effect is small for Q 2 > 3 GeV 2 . 
The average transverse size (p 2 ) g at x ~ 10~ 3 decreases 
by - 15% between Q 2 = 3 and 10 4 GeV 2 , while the effec- 
tive value of a' g in this x-region drops by about half. We 
note that the J/tp electroproduction data from HERA 
[1, Q provide some indication that the effective a' g may 
be smaller than in photoproduction, although the re- 
sults are not fully conclusive. In any case, the a' g in 
the parametrization Eqs. ©-© can be considered as an 
upper limit at values of Q 2 > 10 GeV 2 , as are of interest 
in the applications here. 

Comparatively little is known about the transverse dis- 
tribution of singlet quarks (q + q) at small x. Com- 
parison of the HERA deeply-virtual Compton scatter- 
ing [IU and J/tp production data indicates that singlet 
quarks at x < 10~ 2 are distributed over a larger trans- 
verse area than the gluons, in qualitative agreement with 
theoretical arguments based on the pion cloud contri- 
bution to the parton densities at large b [l9j]. In the 
applications here we are concerned with gluon-induced 
processes; parametrizations similar to Eq. ([5]) could be 
formulated also for the quark distributions. 



III. IMPACT PARAMETER DISTRIBUTION OF 
PROTON-PROTON COLLISIONS 

Using the information on the transverse spatial dis- 
tribution of partons in the nucleon, one can infer the 
distribution of impact parameters in pp collisions with 
hard parton-parton processes. While not directly observ- 
able, the latter determines the spectator interactions in 
such collisions and thus can be studied indirectly through 
measurements of the correlation of hard processes with 
final-state properties. The hard parton-parton process 
is effectively pointlikc in transverse space compared to 
the typical scale of variation of the transverse distribu- 
tions of partons in the colliding hadrons. The impact 
parameter distribution of pp events with a hard gluon- 
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FIG. 3. Overlap integral of the transverse spatial parton dis- 
tributions, defining the impact parameter distribution of pp 
collisions with a hard parton-parton process, Eq. (|10|) . 



gluon process is thus given by the normalized geometric 
probability for two gluons to collide at the same point in 
transverse space: 



P 2 (a;i,X2,&|<3 2 



Pi 



d 2 p2 



Pi + Pi) 



x F g (x uPl \Q 2 ) F g (x 2lP2 \Q 2 ), (10) 

where b = \b\ is the pp impact parameter and p\^, = |pi,2| 
the transverse distances of the two gluons from the center 
of their parent protons (see Fig. [3]) [l|. It satisfies the 
normalization condition J d 2 bP2(xi, x 2 , b\Q 2 ) — 1. With 
the parametrizations of Eq. the convolution integral 
in Eq. (fTO)) can easily be evaluated analytically. In the 
case of symmetric collisions (x = x\ = x%) we find 



PiixMQ 1 



(4ttB 9 )- 1 exp[-b 2 /(4B g )}, 
[m 2 /(127r)] K&/2) 3 iv 3 (m 9 6), 



(11) 



where the parameters B g and m g are taken at the appro- 
priate values of x and Q 2 . 

The impact parameter distribution in minimum-bias 
inelastic pp collisions can be inferred from the pp elastic 
scattering amplitude, which incorporates the information 
on the pp total cross section through the unitarity rela- 
tion. It is given by 

P in (s,b) = [l-|l-r( s ,6)| 2 ]/ CTin ( s ), (12) 

where r(s, b) is the profile function of the pp elastic am- 
plitude in the conventions of Ref. [l| and the inelastic 
cross section <7in(s) is given by the integral J d 2 b of the 
expression in brackets, such that J d 2 b Pj n (s,6) = 1. For 
the purpose of the present study we employ a simple an- 
alytic parametrization of the profile function which satis- 
fies unitarity and reflects the approach to the black-disk 
regime (r — ¥ 1) at small impact parameters: 



T(s,b) = T cxp{-6 2 /[2J3(s)]}, 



(13) 



where ro = 1 and the slope parameter is given in terms 
of the total cross section as B(s) = <7tot(s)/(47r); the in- 
elastic cross section for this profile is <7i n (s) = 3wB(s). 
For the total cross section we use the ex trap olation sug- 
gested by the COMPETE Collaboration [if, which gives 
B(s) = 20.2 (22.8) GeV" 2 at y/s = 7 (14) TeV. The un- 
certainty in the profile function at LHC energies is dom- 
inated by that of the total cross section. The impact 
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FIG. 4. Impact parameter distributions of inelastic pp col- 
lisions at ■yjs = 7 TeV. Solid (dashed) line: Distribution 
of events with a dijet trigger at zero rapidity, 1/1,2 = 0, cf. 
Eq. ((IT), for p T = 100 (10) GeV . Dotted line: Distribution 
of minimum-bias inelastic events, cf. Eq. (|12p . 



FIG. 5. Median impact parameter 6(median) of events with 
a dijet trigger, as a function of the transverse momentum pt, 
cf. Fig. |4] Solid line: Dijet at zero rapidity 2/1,2 = 0. Dashed 
line: Dijet with rapidities 2/1,2 = ±2.5. The arrow indicates 
the median b for minimum-bias inelastic events. 



parameter distributions calculated with Eq. (|13p provide 
a fully satisfactory representation of those obtained with 
more elaborate parametrizations of the pp elastic ampli- 
tude, see Fig. 1 of Ref. [l3[ and references therein. 

Using the above expressions we can now study the in- 
fluence of the trigger conditions on the impact parame- 
ter distribution of pp events at the current LHC energy 
y/s = 7 TeV. The present experiments typically consider 
a jet trigger near zero rapidity, y\ 0, and study the 
characteristics of the underlying events as a function of 
the transverse momentum px of the highest-momentum 
particle in the pseudorapidity interval —2.5 < 77 < 2.5. 
In this setting one integrates over the energy of the bal- 
ancing jet (as well as that of other jets which might arise 
from higher-order processes), which effectively amounts 
to integrating over the momentum fraction of the second 
parton, X2, at fixed x\. Since the distribution is sym- 
metric in the rapidity of the balancing jet, y 2 , and the 
variation of the transverse distribution of partons with x 
is small, cf. Eqs. ([5])-© and Fig. [3J we can to a good ap- 
proximation set ?/2 = and thus take X1.2 at the average 
point 

xi = x 2 = IptI^Ts. (14) 

The scale at which the parton densities are probed is 
of the order Q 2 ~ p T , with a coefficient which remains 
undetermined at leading-order accuracy. Generally, we 
expect the impact parameter distribution in events with 
such a jet trigger to become narrower with increas- 
ing px, because the transverse distribution of partons 
shrinks both with increasing a; 1.2 and with increasing Q 2 . 



The impact parameter distributions with a jet trigger of 
Pt = 10 and 100 GeV are presented in Fig. (U Shown are 
the results obtained with the exponential parametriza- 
tion of Eq. (fTTj) and Eqs. ©-(©J the dipole form leads 
to comparable results. One sees that the change of the 
width of this distribution with px is rather small, because 
the transverse distribution of gluons changes only little 
with x in the range explored here; account of the Q 2 de- 
pendence of the transverse distribution of gluons would 
lead to an additional small change. One also sees that 
the impact parameter distributions with the jet trigger 
are much narrower than that in minimum-bias inelas- 
tic events at the same energy. This quantifies the two- 
scale picture of transverse nucleon structure summarized 
in Fig. [TJ 

The median impact parameter in dijet events, de- 
fined as the value of b for which the integral over P 2 
reaches the value 1/2, is shown in Fig. [5] as a function of 
Pt- For the parametrizations of Eq. pip it is given by 
fr(median) = 1.67 \jB g and 3.08 m" 1 , respectively. The 
results obtained with the exponential and dipole form 
factors differ only by a few percent if the parameters are 
related by Eq. ((3J), indicating that the uncertainty re- 
sulting from our imperfect knowledge of the shape of the 
transverse spatial distribution of gluons is small. The 
uncertainty in fe(median) resulting from the uncertainty 
of B g o in the parametrization Eqs. ©-((!]) is less than 
±10% at pt ~ few GeV. It is seen that the median b 
in jet events drops only very weakly as a function of pt 
for all values above ~ 2 GeV. We estimate that account 
of the Q 2 dependence of the transverse distributions due 
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to DGLAP evolution would change the results in Fig. [5] 
by less than ~ 5%. Also shown is the median b with a 
trigger on a jets at non-zero rapidity y\ = —7/2 = 2.5, 
which amounts to an effective increase of x\ y 2 by a factor 
coshy s» 6, cf. Eq. (fT6)l and the discussion in Sec. [V] In 
all cases, the median impact parameter in jet events is 
far smaller than that in minimum-bias collisions, which 
is given by ^(median) = 1.32 ^/B for the parametrization 
ofEq. 

To conclude this discussion, a comment is in order con- 
cerning the interpretation of the impact parameter dis- 
tributions in pp events with hard processes. Our analysis 
based on Eq. (flOl) shows that pp events with at least one 
hard process (and no other requirements) are on average 
more central than minimum-bias inelastic events. This 
statement concerns the relative distribution of impact pa- 
rameters in a collective of inelastic pp events and how it 
is changed by imposing the requirement of a hard pro- 
cess. One should not confuse this with statements about 
the absolute probability for a hard process (in a certain 
rapidity interval) in a pp collision at certain impact pa- 
rameters. In fact, the analysis of Refs. [U [23] shows 
that there can be a substantial absolute probability for 
a hard process in pp collisions at large b, and that uni- 
tarity places non-trivial restrictions on the dynamics of 
hard interactions in peripheral collisions. 



IV. TRANSVERSE MULTIPLICITY AS AN 
INDICATOR OF HARD DYNAMICS 

The estimates of the previous section show that pp 
events with a hard parton-parton collision are much more 
central than minimum-bias events, and that the average 
impact parameters change only very little for pt above 
~ 2 GeV. At the same time, it is known that the overall 
event characteristics, such as the average multiplicity, de- 
pend strongly on the centrality of the underlying pp col- 
lision. Combining these two observations, we can devise 
a practical method to determine down to which values 
of pt mid-rapidity particle production is predominantly 
due to hard parton-parton collisions. The observable of 
interest is the transverse multiplicity, measured in the 
direction perpendicular to the transverse momentum of 
the trigger particle or jet. It is not directly affected by 
the multiplicity associated with the trigger or balancing 
jets, but is indirectly correlated with the presence of a 
hard process because of its dependence on the centrality. 

Based on the results of Figs.|4]and[5]we predict that the 
transverse multiplicity should be practically independent 
of pt of the trigger as long as the trigger particle orig- 
inates from a hard parton-parton collision which "cen- 
ters" the pp collision. Furthermore, the transverse multi- 
plicity in such events should be significantly higher than 
in minimal-bias inelastic events, since the known mecha- 
nisms of particle production — minijet interactions, mul- 
tiple soft interactions, etc. — are much more effective in 
central collisions. When measuring the transverse multi- 
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FIG. 6. Schematic illustration of the expected dependence of 
the transverse multiplicity, N(pr), on the pr of the trigger. 



plicity as a function of pt of the trigger, we thus expect 
it to increase from its minimum-bias value at low pt and 
become approximately constant at pt ~ few GeV (see 
Fig. [5]). The point where the transition happens, pr.crit, 
indicates the critical value of pt above which particle pro- 
duction is dominated by hard parton-parton processes. 

Interestingly, the predicted increase and eventual flat- 
tening of the transverse multiplicity agrees well with the 
pattern observed in the existing data. At y/s — 0.9 TeV 
the transition occurs approximately at pr.crit ~ 4 GeV 
@, at y/s = 1.8 TeV at p T ,crit ~ 5 GeV 0], and the pre- 
liminary data at 7 TeV indicate somewhat larger values 
of PT.crit = 6 — 8 GeV [1, 0| • We thus conclude that the 
minimum pt for hard particle production increases with 
the collision energy. Note that we consider here an inclu- 
sive trigger; the procedure adopted in the experimental 
analysis (selection of the fastest particle in the measured 
rapidity interval) somewhat enhances the contribution of 
soft mechanisms in particle production. 

It is worth noting that the overall pattern described 
here is reproduced by the tunes of current MC models; 
cf. the comparisons in Refs. 0-0] ■ This is because these 
models effectively include the key feature used in our 
analysis — the narrow impact parameter distribution of 
dijet events (although (b 2 ) in these models is too small by 
a factor ~ 2), and impose a cutoff on the minimal pt of 
the minijets. Our point here is that the observed pattern 
can be explained naturally on the basis of the transverse 
geometry of pp collisions with hard processes, without in- 
volving detailed models. This allows one to determine in 
a model-independent way where the dominant dynamics 
in particle production changes from soft interactions to 
hard parton-parton processes. 

For pt lower than pT,crit the relative contribution of 
hard processes to particle production starts to decrease. 
In terms of the transverse geometry, this means that the 
observed trigger particle can, with some probability, orig- 
inate from either peripheral or central collisions in the 
sense of Fig. [TJ We can estimate the fraction of particles 
produced by hard interactions in this "mixed" region in a 
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simple two-component model, based on the observation 
that the effective impact parameters in soft collisions are 
much larger than those in hard events and do not change 
much with transverse momenta of the produced particles 
[? ]. Thus, we assume that: (i) a trigger particle ob- 
served at given pr originated with a probability Ahard(pT) 
from a hard process, and with probability 1 — Ahard(PT) 
from soft interactions; (ii) the average impact parame- 
ters in both classes of collisions do not depend of the pt 
of the trigger. This allows us to write the pt dependence 
of the transverse multiplicity in the form 

N(p T ) = A hard (p T )7V hard + [1 - A hard (p T )]iV soft , (15) 

where A^ard and N so ft are independent of px- Assum- 
ing that for some sufficiently small px cutoff Ahard(PT) 
is close to zero, we can determine N so {t, which corre- 
sponds to the minimum-bias impact parameter distribu- 
tion presented in FigJH and use it to determine Ahard (pt) 
for pt smaller than pr,crit(s) via Eq. (TT5|) . The data 
indeed indicate that Ahard ^ A^oft; so in the region 
of pr where N (pt) / > 1/3 our estimate is not 

sensitive to the exact value of N so f t . By inspection of 
the data we conclude that the contribution of the hard 
mechanism drops to about half of the total yield for 
p T 1.5 -2,2- 2.5,3 -4GeV for y/s = 0.9, 1.8, 7TeV. 

It is also of interest that for pt > Pt, crit the trans- 
verse multiplicity appears to increase with the collision 
energy faster than the average multiplicity [1, 0] • In the 
leading-twist approximation the perturbative contribu- 
tion is proportional to the product of the gluon densities 
at small x\ t 2 and thus scales as (y/s) 2X , where A is the ex- 
ponent of the gluon density, xg(x, Q 2 ) cx x~ A , and takes 
on values A ~ 0.2 — 0.3 in the Q 2 region of interest here. 
This is roughly consistent with the factor ~ 2 increase 
of the observed transverse multiplicity between y/s = 0.9 
and 7 TeV, which suggests scaling as (y/s) ' 34 [? ]. We 
note that at very small values of x\ or xi the leading- 
twist approximation breaks down because of the onset 
of the black-disk regime in hard interactions, which gen- 
erates a new dynamical scale in the form of the gluon 
density per transverse area; see Ref. |l| for an estimate 
of the relevant values of X\^ and pt- 



V. RAPIDITY DEPENDENCE AS A TEST OF 
UNIVERSALITY 

The basic idea of our approach is that hard processes 
influence the event characteristics by selecting pp colli- 
sions with small impact parameters. Further interesting 
tests of these geometric correlations can be performed by 
measuring the dependence of event characteristics on the 
rapidity of the jets in the trigger. 

In production of dijets at non-zero rapidity only part of 
the center-of-mass energy of the colliding partons is con- 
verted to transverse energy, allowing one to probe larger 
momentum fractions x\.2 at the same pr- For jets with 



symmetric rapidities y\ = — yi = y, 

xi = x 2 = (2p T coshy)/Vs. (16) 

Because partons with larger x±.2 sit at smaller transverse 
distances, the average impact parameters in pp collisions 
with a dijet trigger decrease with increasing y; however, 
the effect is small (see FigE}. Observing the approximate 
y-independence of the transverse multiplicity would test 
that the selection of central collisions does not depend 
on the details of the hard process. Beyond that, we pre- 
dict a small increase in the transverse multiplicity if y 
is increased away from zero at fixed pt- In particular, 
such measurements could separate the effects of the x— 
and Q 2 -dependence of the transverse distribution of par- 
tons on the average impact parameters. At lower pt, the 
dependence of the transverse multiplicity on the rapid- 
ity would help to distinguish between the minijet and 
soft mechanisms of hadron production, as minijets are 
much more centered at small rapidities, while typical soft 
multi-ladder-type interactions lead to correlations over 
large rapidity intervals. An additional advantage of mea- 
surements with the y ^ trigger is that the difference 
between the transverse, forward, and away-side regions 
of particles produced at mid-rapidity is much smaller 
than for the y = trigger. 

The selection of central pp impact parameters by hard 
parton-parton processes could in principle be verified 
not only through the transverse multiplicity, but also by 
measuring event characteristics in rapidity regions which 
are not directly affected by the jet fragmentation of the 
partons in the trigger process. The identification of the 
jet fragmentation regions requires detailed modeling be- 
yond the scope of the present investigation. Assuming 
that one could reliably remove the fragmentation regions, 
several types of interesting correlation measurements be- 
come possible. First, we predict that in the remaining 
rapidity region the multiplicity should become, on aver- 
age, isotropic in the transverse direction, i.e., it should 
attain the value of the previously considered transverse 
multiplicity in all directions and be substantially higher 
than in minimum bias events. Second, this multiplicity 
should not depend on the rapidity of the trigger, y, if 
Pt > PT,crit- Both measurements would directly attest 
to the universality of particle production in central pp 
collisions. 

The present LHC experiments use the central detectors 
to study the underlying event structure in the produc- 
tion of high-pT particles in the pseudorapidity interval 
r\ = ±2.5 (corresponding roughly to the same range in 
rapidity proper, y); the measurements may be extended 
to -5 < 7] < 5 (CMS) and -4.9 < r] < 4.9 (ATLAS) 
using forward detectors. In production of two jets at 
yi = —y2 = y w 2, assuming a rapidity interval of ap- 
proximately ±0.7 for the fragmentation of either jet, the 
rapidity region ±1 should be free of direct jet fragments 
and could be used for the envisaged multiplicity mea- 
surements. Alternatively, one may consider a pair of jets 
at the same positive rapidity, y\ — y 2 = y > 0, and 
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study the multiplicity in the negative rapidity region as 
a function of y. The latter choice would have the advan- 
tage that the parton momentum fractions x\p, change in 
different directions when increasing y from zero, compen- 
sating the effect on the width of the impact parameter 
distribution to first order. 

An interesting phenomenon should occur when extend- 
ing such measurements with symmetric jets at y\ = 
— 2/2 = V to larger rapidities. As discussed in Sec. HH the 
transverse size of the parton distribution decreases with 
increasing x. This leads to a seemingly paradoxical pre- 
diction, that the larger the rapidity interval yi — 2/2 = 2y 
between the jets, the larger the multiplicity in the mid- 
rapidity region. In other words, one expects long-range 
correlations in rapidity which are becoming stronger with 
increase of the rapidity interval. However, the effect is 
rather small; at pt — 5 GeV the median pp impact pa- 
rameter changes from 0.66 fm at y = to 6.1 (0.55) fm at 
y = 2.5 (5). Therefore the ability to measure over a wide 
range of pseudorapidities |^| < 5 would be very helpful 
for studying this effect. 

Deviations from the predicted universality can arise as 
a result of spatial correlations between partons involved 
in the hard collisions and those participating in spectator 
interactions [23[ ■ At the relatively small ^-values probed 
with the central detectors at LHC (x < 10 -2 ), such cor- 
relations are likely to depend weakly on a; and would not 
significantly affect the rapidity dependence. In princi- 
ple, the study of these deviations from universality may 
provide a new window on correlations in the nucleon's 
partonic wave function. 



VI. SUGGESTIONS FOR FUTURE 
MEASUREMENTS 

In addition to the studies described in Sees. II VI and W\ 
several other kinds of measurements could further explore 
the proposed connection between hard processes and the 
transverse multiplicity, or use it to investigate interesting 
aspects of QCD and nucleon structure. 

Energy dependence of transverse multiplicity. Mea- 
surement of the energy dependence of the transverse mul- 
tiplicity in jet events would, in effect, reveal the energy 
dependence of the average multiplicity in central pp colli- 
sions, which is of interest beyond the specific applications 
considered here. In order to avoid change of the average 
impact parameters due to the ^-dependence of the trans- 
verse distribution of partons, it would be desirable to 
compare data at the same values of xi } 2- When increas- 
ing the collision energy from y/s to ^/s one thus needs 
to trigger on jets with rapidities scaled by (1/2) ln(s/so). 

Double dijet trigger. Further reduction of the effective 
impact parameters in pp collisions can be achieved with 
a trigger on multiple hard processes [l|. In the absence 
of transverse correlation between partons, the effective 
impact parameter distribution of events with two dijets 
would be given by the (properly normalized) product of 



distributions P2 in Eq. (fT0| . In the simplest case of two 
dijets with x\ = X2 = X3 — X4, the median b in such 
events would be 1.18 yJB g (1.97 nig 1 ), to be compared 
to 1.67 y/Bg (3.08 m" 1 ) for single jet events. Account 
of transverse correlations between partons reduces the 
difference in effective impact parameters by about half 
[l[. In all, we expect a 15— 20 % reduction of the median b 
with a double dijet trigger, which should manifest itself in 
a further increase of the transverse multiplicity compared 
to single dijet events. 

Other centrality triggers. Knowledge of the depen- 
dence of the transverse multiplicity on b would allow one 
to calibrate other triggers on central pp collisions. In par- 
ticular, it would be interesting explore triggers related to 
particle production in the nucleon fragmentation regions, 
for example leading neutrons. Ultimately one aims here 
at designing a trigger on ultra-central pp collisions, in 
which the effective gluon densities would be comparable 
to those reached in heavy-ion collisions [24j . 

Quark vs. gluon-induced jets. It would be interesting 
to compare the transverse multiplicities and other under- 
lying event characteristics for quark-antiquark induced 
hard processes like W ± , Z° production and gluon-gluon 
induced processes. Another possibility would be to con- 
sider large— 1?7| dijet production and separate quark- and 
gluon induced jets using the different jet shapes. This 
would allow one to observe a possible difference between 
the transverse distributions of quarks and gluons, com- 
plementing studies of hard exclusive processes in ep/jp 
scattering. 

Deconvolution of impact parameter dependence. If the 
fluctuations of the multiplicity between events at a given 
impact parameter are not too large, one can attempt a 
program of deconvolution of the impact parameter de- 
pendence of the multiplicity, using information on the 
impact parameter dependence of dijet and 4-jet rates 

VII. SUMMARY AND DISCUSSION 

The transverse multiplicity in pp collisions with a jet 
trigger provides an effective way to determine the rate 
of processes initiated by hard parton-parton interactions 
down to rather small transverse momenta. Further anal- 
ysis of the data can then establish whether the observed 
particle production rates are consistent with perturba- 
tive QCD predictions. Studies of the dependence of the 
transverse multiplicity on the transverse momenta and 
rapidities of the jets, and on the collision energy, can 
provide a better understanding of the impact parameter 
dependence of the underlying event characteristics and 
allow one to refine the use of the transverse multiplicity 
as an indicator of the dynamics in particle production. 

The results of the present study have implications also 
for the use of hard processes in new particle searches. 
In the production of a Higgs boson with a mass Mh ~ 
100 GeV the parton momentum fractions X1.2 are the 
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same as in dijet events at zero rapidity with pt = 
Ma/ 2 = 50GcV. The two types of events thus involve 
the same average pp impact parameters. This allows 
one to describe the background to new particle produc- 
tion processes much more accurately than on the basis of 
minimum-bias event characteristics. 

Our analysis relies crucially on information about 
the transverse spatial distribution of gluons from exclu- 
sive J /ip photo/electroproduction and similar processes. 
While the region x < 0.01 has been covered by HERA, no 
data of comparable precision are available at larger x (see 
Sec. |TTJ) . This is unfortunate, as the production of new 
particle with masses ~ 10 2 GeV requires partonic colli- 
sions at precisely such momentum fractions. The region 
x > 0.01 is also particularly interesting for nucleon struc- 
ture studies [13, [3 ■ The J/tp results expected from the 
COMPASS experiment, as well as measurements with a 
future Electron-Ion Collider (EIC), thus would have a 
major impact on both areas of study. 
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